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NOTE  ON  COMPUTER  GRAPHICS  FOR 
MAXIMUM  ENTROPY  SPECTRAL  ANALYSIS 


Summary i  Computer  graphics  provides  an  effective  three-dimensional 
display  of  power  spectral  density  for  varying  parameters.  Based  on 
the  software  we  have  developed,  the  computer-graphics  display  is 


presented  for  Burg’s  maximum  entropy  spectral  analysis  and  Fougere's 
(nonlinear)  maximum  entropy  spectral  analysis.  The  significantly 
greater  dynamic  range  of  the  Fougere's  spectrum  clearly  illustrates 
the  superiority  of  the  Fougere’s  method.  The  graphics  software 
listing  is  provided  in  the  Appendix  of  the  report. 

Spectral  Display  of  Sunspot  Numbers 

Fig.  la  is  a  plot  of  January  Zurich  sunspot  numbers  for  200 
data  points  for  the  period  1779  to  1978.  Fig.  lb  is  the  display  of 
Burg's  spectrum  of  Fig.  la  for  filter  weights  10  to  58  and  frequency 
interval  0.05  to  0.15  cycles  per  year.  All  spectrum  plots  are  in 
logarithmic  scale  with  base  10.  The  three  numbers  on  the  top  of 
the  photo  are  A=  minimum  spectral  magnitude,  B=  maximum  spectral 
magnitude,  and  C=  difference  between  B  and  A,  in  log  scale. 

For  Fig.  lb,  A=2.798,  B=5*113.  C=2.315.  In  all  graphical  displays 
of  the  report,  450  points  are  used  for  computation  in  frequency  and 
weight  number  axes  even  though  only  225  points  are  actually  used 

for  display  (See  Appendix). 

Fig.  2a  is  a  section  of  25  data  points  from  Fig.  la  to  cover 
the  period  1954-1978.  Fig.  2a  is  the  Burg's  spectra  for  filter 
weights  of  1  to  17.  The  dynamic  range  is  A=6.28l,  B=l4.424,^ 
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C=8.l43.  Fig.  2c  is  the  Fougere's  spectra  of  Fig.  2a  for  filter 
weights  of  1  to  17.  The  dynamic  range  is  A =4. 403,  B=l?.225  and 
C=12.822.  The  increase  in  dynamic  range  by  using  the  Fougere's 
method  is  about  46,8dB. 

Spectral  Display  of  Sinewaves 

Figure  3a  is  the  sinewave  considered  which  is  given  by  the 
equation 

x(t)=  sin (2  ft  +0)  +  SNR  x  n(t) 

with  DT( sampling  period )=0. 05 »  f=l,  0=45°,  SNR=50dB,  and  n(t)  is 
zero  mean, unit  variance  Gaussian  noise.  25  data  points  are  used 
in  spectral  analysis.  Fig.  Jo  is  the  Burg's  spectra  for  filter 
weights  of  1  to  17  and  frequency  interval  of  0  to  2.7  Hz.  The 
dynamic  range  is  A=-l4.l08,  B=3.085,  C=17.193»  Fig.  3c  is  the 
Fougere's  spectra  for  the  same  filter  weights  and  frequency 
interval.  The  cynamic  range  is  A=-l4.584,  B=10.353,  C=24.93 7. 

The  increase  in  dynamic  range  of  the  spectral  magnitude  by 
using  the  Fougere's  method  is  about  77»5dB. 

Useful  References: 

1.  C.  H.  Chen,  "Spectral  resolution  of  Fougere's  maximum  entropy 
spectral  analysis",  Proc.  of  the  IEEE,  June  1981. 

2.  C.  H.  Chen,  J.  Chen  and  C.  Yen,  "A  Minicomputer  implementation 
of  Fougere's  maximum  entropy  spectrum  analysis  method", 
Technical  Report,  SMU-EE-TR-80-7,  August  20,  1981. 
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